Re: dependency diagram for druckprag

[email protected] Sat, 25 Feb 2012 20:35:10 -0600
Newsgroups gmane.comp.mathematics.tochnog.user
Message-ID <[email protected]>
--Apple-Mail=_73B801B6-E243-427A-8738-10C74240DB18
Content-Transfer-Encoding: quoted-printable
Content-Type: text/plain;
	charset=windows-1252

Ning,

I read your comment and still believe that you are not entering the =
dependency diagram items the proper way.=20
I checked the code and the open source version allows to enter more than =
one property. See the attached note=20
from the manual on dependency_item:

"=85.
Temperature dependent Young's moduli in two maxwell chains of =
element_group 1 ( for the first chain the=20
Moduli 1.e10, 1.e9,...  for the second chain the moduli 1.e12, 1.e11,... =
 all relaxation times are 1.e-2 . ):

dependency_item 1 -group_materi_maxwell_chain 1 -temp 4=20
dependency_diagram 1=20
1. 2. 3. 4.=20
1.e10 1.e9 1.e8 3.e5=20
1.e-2 1.e-2 1.e-2 1.e-2=20
1.e12 1.e11 1.e10 3.e7=20
1.e-2 1.e-2 1.e-2 1.e-2=20
=85"

As you can see, you enter the list of unknowns, in the example above, =
the values of -temp, then you enter the list=20
of values of the
first chain moduli, and then the list of the  relaxation time values, =
then the second list of the second chain moduli,
and then the list of the second relaxation times for the secon chain.

for your case, upon reviewing your file, you have=20

"=85
group_materi_plasti_druckprag 1 0.611 1.5e6 0.175
group_integration_points 1 -maximal

dependency_item 1 -group_materi_plasti_druckprag 1 -kap 2
dependency_diagram 1 0. 0.1 0.611 1.5e6 0.175 0.611 5246.89 0.175
=85"

Therefore, I understand that  you are using:
group_materi_plasti_druckprag index  phi c phiflow

and that your inital your initial values are phi=3D 0.611, c =3D 1.5e6 =
and phiflow=3D 0.175

You are keeping phi and phiflow constant and are modifying c. Your =
dependency diagram should therefore be as follows:

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

dependency_item 1 -group_materi_plasti_druckprag 1 -kap 2
dependency_diagram 1 0. 0.1  0.611 0.611 1.5e6 5246.89 0.175 0.175

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

So I am giving two values of -kap, two values of phi, two values of c, =
and two values of phiflow.

I changed the values and it runs well and very fast. I am attaching a =
copy of the modified file for your review.

Just a couple of reminders:=20

1.	The updated, that is current version of tochnog allows the use =
of the Lapack solver which is faster than the=20
	default biconjugate solver. I built the version for Mac OSX but =
currently do not have access to Linux or Windows compilers.

2.	Rather than using "group_integration_points 1 -maximal" you can =
use the modern version of tochnog: "number_of_integration_points N" ,
	where N is any number of integration points that you want to =
use. If you read the manual you will learn that  -maximal has a limit,
	 but "number_of_integration_points" does not. Just be reasonable =
and understand how many points you will need.
3.	If you use "control_print_gid "number"  -yes"  and use the same =
number or index for n as for control_timestep, then you can print=20
	and plot (at least in GiD) the results for any time step and not =
only at the end of the calculations.
4.	Finally, I keep suggesting to use the automatic iterations =
option in tochnog, as it allows you to control the reslution or error  =
of the=20
	problem. I included that in the modified file attached.

Good luck and please let me know if this works.

Best regards=

--Apple-Mail=_73B801B6-E243-427A-8738-10C74240DB18
Content-Disposition: attachment;
	filename=Drucker.dat
Content-Type: application/octet-stream;
	name="Drucker.dat"
Content-Transfer-Encoding: quoted-printable

echo=20-yes=0Anumber_of_space_dimensions=202=0Amateri_velocity=0A=
materi_stress=0Amateri_strain_plasti=0Amateri_plasti_kappa=0A=
number_of_integration_points=204=0Aend_initia=0A=0Anode=201=200.000=20=
0.000=0Anode=202=200.500=200.000=0Anode=203=201.000=200.000=0Anode=204=20=
0.000=200.500=0Anode=205=200.500=200.500=0Anode=206=201.000=200.500=0A=
node=207=200.000=201.000=0Anode=208=200.500=201.000=0Anode=209=201.000=20=
1.000=0A=0Aelement=201=20-quad4=201=202=204=205=0Aelement=202=20-quad4=20=
2=203=205=206=0Aelement=203=20-quad4=204=205=207=208=0Aelement=204=20=
-quad4=205=206=208=209=0A=0Aelement_group=201=201=0Aelement_group=202=20=
1=0Aelement_group=203=201=0Aelement_group=204=201=0A=0Astart_define=0A=20=
=20bottom_edge=20geometry_line=201=0Aend_define=0A=0Astart_define=0A=20=20=
left_edge=20geometry_line=202=0Aend_define=0A=0Astart_define=0A=20=20=
top_edge=20geometry_line=203=0Aend_define=0A=0Astart_define=0A=20=20=
right_edge=20geometry_line=204=0Aend_define=0A=0Abottom_edge=200.0=200.0=20=
1.0=200.0=200.000001=0Aleft_edge=200.0=200.0=200.0=201.0=200.000001=0A=
top_edge=200.0=201.0=201.0=201.0=200.000001=0Aright_edge=201.0=200.0=20=
1.0=201.0=200.000001=0A=0Abounda_unknown=200=20-bottom_edge=20-vely=0A=
bounda_unknown=201=20-left_edge=20-velx=0Abounda_unknown=202=20-top_edge=20=
-vely=0Abounda_time=202=200.0=20-0.01=201.e6=20-0.01=0A=0Agroup_type=201=20=
-materi=0Agroup_materi_elasti_young=201=2025e9=0A=
group_materi_elasti_poisson=201=200.3=0Agroup_materi_memory=201=20=
-updated_without_rotation=0Agroup_materi_plasti_druckprag=201=200.611=20=
1.5e6=200.175=0A=0A(options_solver=20=20-matrix_lapack)=0A=0A=
dependency_item=201=20-group_materi_plasti_druckprag=201=20-kap=202=0A=
dependency_diagram=201=200.=200.1=20=200.611=200.611=201.5e6=205246.89=20=
0.175=200.175=0A=0Acontrol_timestep=2010=200.01=205.5=0A=
control_timestep_iterations=20=20=20=20=20=20=20=20=20=20=2010=20=205=0A=
control_timestep_iterations_automatic=2010=200.001=200.01=0D=0A=
control_print=20=20=20=20=2010=20=20-time_current=20=
-post_node_rhside_ratio=0Acontrol_print_gid=2010=20-yes=0A=0Aend_data=0A=

--Apple-Mail=_73B801B6-E243-427A-8738-10C74240DB18
Content-Transfer-Encoding: 7bit
Content-Type: text/plain;
	charset=us-ascii



Fernando Lorenzo
[email protected]




--Apple-Mail=_73B801B6-E243-427A-8738-10C74240DB18
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

------------------------------------------------------------------------------
Virtualization & Cloud Management Using Capacity Planning
Cloud computing makes use of virtualization - but cloud computing 
also focuses on allowing computing to be delivered as a service.
http://www.accelacomm.com/jaw/sfnl/114/51521223/
--Apple-Mail=_73B801B6-E243-427A-8738-10C74240DB18
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

_______________________________________________
Tochnog-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/tochnog-users

--Apple-Mail=_73B801B6-E243-427A-8738-10C74240DB18--