|November 4, 2010, 11:36||
Orthotropic solid displacement
I'm trying to write a solver, loosely based on solidDisplacementFoam, which will work on orthotropic solids. The stiffness constants and the density are specified as fields (and there is also a misfit field f, since what I'm actually trying to do involves coherent interfaces between materials with different lattice constants).
Since the system is anisotropic I can't just calculate Lamé parameters, but since OpenFOAM can't handle the full fourth-rank 3x3 stiffness tensor I can't just write sigma = c : epsilon (where sigma is the stress tensor and epsilon is the strain tensor). However, since it's orthotropic I can work with vectors for the pressure and shear terms for both stress and strain, and two 3x3 tensors with the elastic constants in (for c_ij where i or j are 1, 2 or 3; and for c_kk where k is 4, 5 or 6).
rho * fvm::d2dt2(D) == fvc::div(sigma)
epsilon = 0.5*(fvc::grad(D) + fvc::grad(D)().T());
sigmaPressure = cPressure & epsilonPressure;
sigmaShear = cShear & epsilonShear;
I'm calculating the strain tensor from the displacement, setting up two strain vectors, getting two stress vectors from those via the stiffness tensors, and rebuilding the stress tensor which I'd like to put back on the right-hand side of the equation. But is it even possible to solve a system of equations like this, rather than having the spatial derivatives of D explicitly on the right-hand side?
I appreciate any feedback at all...
|heteroepitaxy, semiconductor, soliddisplacementfoam, strain|
|Thread||Thread Starter||Forum||Replies||Last Post|
|How model radiation of both solid wall and fluid£¿||Harry Qiu||FLUENT||2||February 3, 2013 23:04|
|how to use Total Mesh Displacement for making a subject translating and rotating?||snailstb||CFX||3||September 16, 2009 04:06|
|Defining Solid Zones Not working||elmcmaster||FLUENT||3||May 22, 2009 07:14|
|Multiple Solid Domains - Interfaces||Scott||CFX||8||July 31, 2008 15:20|
|CFX4.3 -build analysis form||Chie Min||CFX||5||July 12, 2001 23:19|