MathGroup Archive 2002

[Date Index] [Thread Index] [Author Index]

Search the Archive

RE: Solve of inverse quadratic finite element problem - resend

  • To: mathgroup at smc.vnet.net
  • Subject: [mg36481] RE: [mg36453] Solve of inverse quadratic finite element problem - resend
  • From: "DrBob" <drbob at bigfoot.com>
  • Date: Sun, 8 Sep 2002 03:31:20 -0400 (EDT)
  • Reply-to: <drbob at bigfoot.com>
  • Sender: owner-wri-mathgroup at wolfram.com

If x1 through x20 aren't unknowns, you should have told me their values.
If you can't, they are UNKNOWN.

True, you're not solving for them... but they affect the dimensionality
of the problem just as if you WERE solving for them.  You have a
20-dimensional space full of contingencies -- solution forms that depend
on the values of x1 through x20.  The Solve function won't deal with
contingencies even if it could, and if it tried, there would be too
many.

Bobby

-----Original Message-----
From: purcell [mailto:chris.purcell at drdc-rddc.gc.ca] 
To: mathgroup at smc.vnet.net
Subject: [mg36481] RE: [mg36453] Solve of inverse quadratic finite element problem
- resend


The constants x1 through x20 are constants.
That leaves 3 equations in 3 unknowns.
Thanks for looking at this. Chris.



At 06:40 PM 9/7/2002 -0500, you wrote:
>You have three nonlinear (fourth-order) equations and 23 unknowns.  A
>faster computer won't help any.
>
>Bobby
>
>-----Original Message-----
>From: purcell [mailto:chris.purcell at drdc-rddc.gc.ca]
To: mathgroup at smc.vnet.net
>Sent: Saturday, September 07, 2002 1:54 AM
>Subject: [mg36481] [mg36453] Solve of inverse quadratic finite element problem -
>resend
>
>
>Would someone with a very fast machine and lots of memory be willing to
>try
>this Solve for me?
>It is the inverse of the 20 node quadratic hexahedral mapping used in
>finite element analysis.
>None of my computers can handle this - they run out of memory (using
>Version 4.2).
>
>(*Hex20 Node definition in global coordinates *)
>Clear[
>x1, y1, z1,
>x2, y2, z2,
>x3, y3, z3,
>x4, y4, z4,
>x5, y5, z5,
>x6, y6, z6,
>x7, y7, z7,
>x8, y8, z8,
>x9, y9, z9,
>x10, y10, z10,
>x11, y11, z11,
>x12, y12, z12,
>x13, y13, z13,
>x14, y14, z14,
>x15, y15, z15,
>x16, y16, z16,
>x17, y17, z17,
>x18, y18, z18,
>x19, y19, z19,
>x20, y20, z20];
>
>(* local coordinates *)
>Clear[u, v, w];
>
>(* Global co-ordinates *)
>Clear[x, y, z];
>
>(* corner nodes *)
>N1= (1-u)*(1-v)*(1-w)*(-2-u-v-w)/8;
>N3= (1+u)*(1-v)*(1-w)*(-2+u-v-w)/8;
>N5= (1+u)*(1+v)*(1-w)*(-2+u+v-w)/8;
>N7= (1-u)*(1+v)*(1-w)*(-2-u+v-w)/8;
>N13=(1-u)*(1-v)*(1+w)*(-2-u-v+w)/8;
>N15=(1+u)*(1-v)*(1+w)*(-2+u-v+w)/8;
>N17=(1+u)*(1+v)*(1+w)*(-2+u+v+w)/8;
>N19=(1-u)*(1+v)*(1+w)*(-2-u+v+w)/8;
>(*  to u nodes *)
>N2= (1-u^2)*(1-v)*(1-w)/4;
>N6= (1-u^2)*(1+v)*(1-w)/4;
>N14=(1-u^2)*(1-v)*(1+w)/4;
>N18=(1-u^2)*(1+v)*(1+w)/4;
>(*  to v nodes *)
>N4= (1+u)*(1-v^2)*(1-w)/4;
>N8= (1-u)*(1-v^2)*(1-w)/4;
>N16=(1+u)*(1-v^2)*(1+w)/4;
>N20=(1-u)*(1-v^2)*(1+w)/4;
>(*  to w nodes *)
>N9= (1-u)*(1-v)*(1-w^2)/4;
>N10=(1+u)*(1-v)*(1-w^2)/4;
>N11=(1+u)*(1+v)*(1-w^2)/4;
>N12=(1-u)*(1-v)*(1-w^2)/4;
>
>(* solve the inverse transform *)
>Solve[{
>x1*N1+x2*N2+x3*N3+x4*N4+x5*N5+x6*N6+x7*N7+x8*N8+x9*N9+x10*N10+
>x11*N11+x12*N12+x13*N13+x14*N14+x15*N15+x16*N16+x17*N17+x18*N18+x19*N19
+
>x20*N20-x==0,
>y1*N1+y2*N2+y3*N3+y4*N4+y5*N5+y6*N6+y7*N7+y8*N8+y9*N9+y10*N10+
>y11*N11+y12*N12+y13*N13+y14*N14+y15*N15+y16*N16+y17*N17+y18*N18+y19*N19
+
>y20*N20-y==0,
>z1*N1+z2*N2+z3*N3+z4*N4+z5*N5+z6*N6+z7*N7+z8*N8+z9*N9+z10*N10+
>z11*N11+z12*N12+z13*N13+z14*N14+z15*N15+z16*N16+z17*N17+z18*N18+z19*N19
+
>z20*N20-z==0},
>{u,v,w}]
>Christopher J. Purcell
>Defence R&D Canada - Atlantic
>9 Grove St., PO Box 1012
>Dartmouth NS Canada B2Y 3Z7
>Tel 902-426-3100 x389, Fax 902-426-9654
>E-mail:    chris.purcell at drdc-rddc.gc.ca

Christopher J. Purcell
Defence R&D Canada - Atlantic
9 Grove St., PO Box 1012
Dartmouth NS Canada B2Y 3Z7
Tel 902-426-3100 x389, Fax 902-426-9654
E-mail:    chris.purcell at drdc-rddc.gc.ca





  • Prev by Date: Re: ListIntegrate Info.
  • Next by Date: RE: Re: fourier transform time
  • Previous by thread: Re: aligning equations
  • Next by thread: RE: Solve of inverse quadratic finite element problem - resend